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We explore the phenomenology of a two Higgs doublet model where both Higgs doublets couple to 
up-type and down-type fermions with couplings determined by the minimal flavor violation ansatz. 
This "2HDM Type MFV" generalizes 2HDM Types I-IV, where the decay rates of /i -> 6b and 
h — >■ T^T~ are governed by MFV couplings independent of the Higgs couplings to gauge bosons 
or the top quark. To determine the implications of the present Higgs data on the model, we 
have performed global fits to all relevant data. Several surprisingly large effects on the light Higgs 
phenomenology can arise: (1) The modified couplings of the Higgs to fermions can enhance the 
/i — > 77 rate significantly in both VBF production (up to a factor of 3 or more) and the inclusive 
rate (up to a factor of 1.5 or more). (2) In the 2HDM Type MFV, the constraints on a light charged 
Higgs are milder than in 2HDM Types I-IV. Thus, there can be substantial charged Higgs loop 
contribution to the di-photon rate, allowing further enhancements of the di-photon rates. (3) The 
h — > r"'"r~ rate can be (highly) suppressed, independently of the other decay channels. Furthermore, 
we studied the correlation between the light Higgs and the heavy Higgs phenomenology. We showed 
that even small deviations from the decoupling limit would imply good prospects for the detection 
of the heavy Higgs boson. In some regions of parameter space, a substantial range of Mh is already 
either ruled out or on the edge of detection. Finally we investigated the possibility that the heavy 
Higgs is close in mass to the light Higgs, providing additional h/H bb rate, as well as confounding 
the extraction of properties of the Higgs bosons. 



I. INTRODUCTION 

How many scalar Higgs doublets are in Nature? The 
intriguing possibility of additional Higgs doublets has 
taken on heightened importance in light of the obser- 
vations by ATLAS [1] and CMS [2] of a particle consis- 
tent with the Higgs boson. One of the simplest exten- 
sions of the Higgs sector of the standard model (SM) is 
a two-Higgs doublet model (2HDM). Glashow and Wein- 
berg long ago recognized that a general 2HDM, in which 
both Higgs doublets couple arbitrarily to the quarks and 
leptons of the standard model, would induce excessively 
large flavor-changing neutral currents (FCNCs) [3] . This 
has led to the various well-known (and less-known) classi- 
fications or "Types" of 2HDM models with natural flavor 
conservation, where the different types of fermions (up- 
type quarks, down- type quarks and leptons) couple to 
one Higgs doublet only (for a recent review see [4]). We 
will refer to these models as 2HDM Types TIV, follow- 
ing [5]. 

In the absence of enlarged symmetries, there is noth- 
ing to enforce exact flavor-conservation of the Higgs cou- 
plings. Thus, it behooves us to understand the full extent 
of what is possible within a general 2HDM, when both 
Higgs doublets couple to all quarks and leptons, but with 
controlled couplings such that FCNC constraints are sat- 
isfied. Various mechanisms have been discussed in the 
literature to protect 2HDMs from too large FCNCs [6- 
16]. Among them, minimal fiavor violation (MFV) [6-9] 
provides a simple ansatz to allow general Higgs couplings 
without excessive FCNCs. 

In this paper we consider a 2HDM in which both Higgs 



doublets couple to quarks and leptons with MFV cou- 
plings. This model was proposed in [6] and it has been 
studied in [17-23] , where the main focus was on the flavor 
phenomenology of the model. In this work we study the 
model in view of the current results from Higgs searches. 
As we will see, the 2HDM Type MFV generalizes the 
various 2HDM Types I-IV that have been recently con- 
sidered in the post-125 GeV literature [5, 24-31]. The 
principle differences between the various Types of 2HDM 
are the couplings to the b quark and the r lepton. In 
the 2HDM Type MFV, these couplings can be modified 
independently from each other, and also independently 
from the Higgs couplings to the top quark and the weak 
gauge bosons. Those aspects of Higgs physics dependent 
on the couplings to the top quark or gauge bosons are 
largely equivalent to the other Types of 2HDM. Hence, 
the production cross section through gluon fusion, vector 
boson fusion (VBF) production, and associated produc- 
tion can be fairly accurately described by 2HDM Types 
TIV couplings (determined by tan/3 and a). By con- 
trast, the 2HDM Type MFV can have completely dif- 
ferent widths Fjf, and Tt-t as compared with any of the 
fiavor-preserving Types of 2HDMs. This has very signif- 
icant effects on Higgs physics, not only on the branching 
fractions h ^ bb and h — > r+r^, but to all modes since 
the total width of the light Higgs boson is dominated by 

There is one additional intriguing possibility, in which 
both light and heavy CP-even Higgs bosons are light (for 
a recent discussion, see [32-36]). The presence of two 
light CP-even Higgs bosons with modified couplings to 
bb and r+r" can potentially allow for increased branch- 



ing fractions to multiple interesting modes. For instance, 
the 77 rate can be substantially larger through one light 
Higgs boson without sacrificing a significant rate into bb 
due to the presence of the second light Higgs boson. The 
hints for rates into bb by the CDF and DO collabora- 
tions [37] could, if this scenario is right, suggest they 
have observed the second Higgs boson. 



II. HIGGS COUPLINGS TO FERMIONS 

The most general couplings of two Higgs doublets Hi 
and H2, with hypercharge 1/2 and —1/2 respectively, to 
the SM fermions have the form 



^ 3 {yu)ij H2QiUj + {vu)i] 

+iye)ij HiL.Ej -(- {yi),j hIL.Ej + h.c 



(1) 



where both Higgs bosons couple to up-type and down- 
type fermions as well as leptons. The resulting masses of 
the fermions are given by 

VSR I 1 - \ VCR I _ \ 

»™« = \Vu + — 2/u J , mrf = — (^yd -I- tfiyd J , 
= ^ [vt + tfiVi) , (2) 

where Vi = vsj^ and vi — vcp are the SU{2)l breaking 
vacuum expectation values (vevs) of the two doublets, 
tan/3 = = V2/V1 is their ratio, and v ~ 246 GeV. 

In general, one cannot take generic y and y couplings 
without inducing huge tree-level contributions to FCNCs. 
One way to avoid tree-level FCNCs in a 2HDM is to im- 
pose an approximate discrete or continuous symmetry to 
the model, where Yukawa couplings to only one Higgs 
doublet are allowed at tree-level. Then, small couplings 
to the other Higgs doublet can arise once loop effects or 
higher dimensional operators are considered [18].^ These 
already small couplings can induce sizable contributions 
to FCNCs. We ensure that FCNCs are avoided by tak- 
ing the couplings y to obey the MFV ansatz [6-9] . In the 
quark sector, the MFV assumption states that there are 
only two spurious that break the global SU{3)^ quark fla- 
vor symmetry of the standard model gauge sector. This 
implies that the couplings y and y are not independent 
of each other, allowing us to write [6] 
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(4) 



with parameters Ci that can in general be complex. The 
terms with e" as well as the other higher order terms 



^ Well known examples are Minimal Supersymmetric Standard 
Model (MSSM) where the jji are loop induced [38-41], or be- 
yond the MSSM (BMSSM) models, where they can arise from 
higher dimensional operators [23, 42, 43]. 



containing both y^ and yd couplings lead to flavor chang- 
ing neutral Higgs couplings, that are nevertheless con- 
trolled by CKM matrix elements and thus naturally 
small. These terms can lead to interesting effects in B 
physics as discussed in [18-20, 22, 23]. Here we are in- 
terested in the impact of the flavor conserving part of 
the Higgs couplings on Higgs collider phenomenology. 
We remark that the higher order terms in the expan- 
sions in Eqs. (3-4) can induce non-universalities between 
the couplings to the first two and the third generation 
of fermions. Since Higgs phenomenology is dominated 
by the third generation, we will set all the higher order 
terms to zero. (In fact, even if the higher order terms 
were generated radiatively [44] , they do not impact Higgs 
phenomenology.) We consider only 



2/u = f-uyu , yd = <^dyd 



(5) 



This is the "aligned 2HDM" framework presented in [12]. 
In the following we concentrate on the CP conserving 
case with real. For the lepton sector, we analogously 
assume that 



yi = £iyi 



(6) 



with a real proportionality factor e^, leading to flavor- 
conserving Higgs-lepton couplings (by contrast, see [45, 
46] for studies of lepton flavor- violating Higgs decays). 
The Ci are flavor-universal, so we can interchangeably 
use the subscripts u ^ t, d ^ b and £ o r in the 
following. 

In the setup outlined above, we can directly express 
the Yukawa couplings in terms of the measured quark 
masses 
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(7) 



In general, the three real parameters are not all phys- 
ical. One of them can always be reabsorbed by a re- 
definition of the two original Higgs doublets. For ex- 
ample, we can choose iJ„ cx H2 + ^uH\ to be the com- 
bination of Higgs fields that couples to up-type quarks. 
This choice of basis for the Higgs doublets corresponds 
to setting e„ — 0, which we assume, without loss of gen- 
erality, throughout the paper. This implies the 2HDM 
Type MFV couplings to the top quark are identical to 
the other 2HDM Types I-IV. Furthermore, this choice of 
basis uniquely defines tan /3 as the ratio of the vev of the 
Hu field with coupling e„ = and the vev of the orthog- 
onal Higgs field Hd (for a discussion of basis invariant 
quantities in 2HDMs see [47-49]). 

We now write the couplings of the Higgs boson mass 
eigenstates with the SM quarks and leptons, as well as 
the gauge bosons. The two Higgs doublets comprise 8 
real scalar fields, three of which are the usual Goldstone 
bosons G and that provide the longitudinal compo- 
nents of the Z and bosons. The remaining physical 
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Higgs bosons consist of two CP-even scalars h and H, 
one CP-odd scalar A and the charged Higgs H^. Their 
interaction Lagrangian is 
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where Vij are elements of the CKM matrix and cw is 
the cosine of the Weinberg angle. For our purposes it is 
justified to take the neutrinos to be massless. For the 
reduced couplings ^, that parameterize the deviations of 
the Higgs couplings from the Yukawa couplings in the 
SM, one finds 
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The angle a diagonalizes the mass matrix of the two CP- 
even Higgs bosons. The couplings of the Higgs bosons 
to the top quark and the gauge bosons are identical to 
the 2HDM Types I-IV. The couplings to the down-type 
quarks and leptons are in general different, parameterized 
by and e^. These couplings interpolate continuously 
between the couplings of the well-studied 2HDM Types I- 
IV, recovering the different Types in the following limits: 
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(10) 



The Higgs couplings satisfy the following sum rules 

i+(e)' = (a' + (e)' 

= 1+4, (11) 

I' a 
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(13) 
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We complete this section with a comment on the free 
parameters of the model and on the Higgs potential: Cru- 
cial parameters for our analysis are and ££, the an- 
gles a and /3 and also the masses of the physical Higgs 
bosons Mil: Mh, Ma, and Mfj±. The masses and an- 
gles are determined by the parameters of the Higgs po- 
tential. As presented in Appendix A, we consider the 
most general 2HDM scalar potential, including also quar- 
tic couplings that are often not considered in the liter- 
ature, e.g., {H^H^Y, {H2Hi)hIHi, and {H2Hi)hIH2, 
c.f. Eq. (Al). Indeed, any symmetry imposed to for- 
bid all of these operators would also forbid the fermionic 
couplings proportional to e to the other Higgs doublet. 
Given the most general structure of the potential allows 
us to treat the Higgs masses Mh, Mh, Ma, and Mjj±, 
as well as the angles a and /3, as free parameters in our 
numerical analysis of Sees. IV and V. This approach is 
justified so long as we are not in the decoupling regime 
M\ ^ w^, where the Higgs masses and the angles are 
strongly correlated: 



Ml 



a 



il4± + 0{X,v^) , (15) 
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(16) 



We explicitly checked that the scenarios discussed be- 
low can be realized by appropriate choices of the quartic 
couplings in the Higgs potential that are compatible with 
constraints from perturbativity and vacuum stability. 



III. HIGGS PRODUCTION CROSS SECTIONS 
AND BRANCHING RATIOS 



For the production and decay of the light Higgs boson, 
there are several changes with respect to the standard 
model rates. Some of these changes do not depend on the 
Type of the 2HDM. This is because, as we saw from the 
previous section, the modification of the top couplings 
as well as the modifications of the gauge boson couplings 
are the same with respect to the various Types of 
2HDMs. This implies several simplifications when com- 
paring one Type of 2HDM with another. Here we wish 
to present the cross sections in the 2HDM Type MFV, as 
well as to elucidate the similarities or differences between 
the different Types of 2HDMs. 
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A. Production Cross Sections 

At the LHC and the Tevatron, the dominant produc- 
tion mechanism for a SM-hke Higgs is gluon fusion. The 
modifications of the couphngs of the Higgs to top and 
bottom quarks result in a modification of this produc- 
tion cross section, that proceeds through top and bot- 
tom loops. In our numerical analysis we use HIGLU [50] 
to compute the various parts of the gluon fusion cross sec- 
tion including top and bottom quarks. Numerically, the 
bottom quark loop contribution is generally quite small. 
In the SM, it contributes at the level of ^-^10% for a Higgs 
boson with = 125 GeV. This also follows for 2HDM 
models because, as we will see in the next section, the 
ATLAS and CMS data do not favor a large enhance- 
ment of the Higgs-bottom coupling. If we consider only 
the top quark contribution, then using the results from 
the previous section, the modification of the top cou- 
pling = Caj sp, is the same among the various Types 
of 2HDMs. Hence, the result for the gluon fusion cross 
section takes a simple (approximate) form 



MFV 

ggh 



2HDM 

ggh 



"ggh.tt 



h\2 



(17) 



where cr^J^ is the cross section in the the 2HDM Type 
MFV, cr^Hp^ is the cross section in the the 2HDM Types 
I- IV, and (Jg^ tt is the top loop contribution of the stan- 
dard model. 

Higgs production through vector boson fusion or in 
association with W or Z bosons plays an important role 
in Higgs searches. In 2HDMs, they scale to an excellent 
approximation with the coupling of the Higgs boson to 
weak vector bosons 
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and thus 
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In our numerical results, we take the SM cross sections 
from the LHC Higgs cross section working group [51, 52]. 

Finally, Higgs boson production in association with top 
or bottom quarks is strongly suppressed in the SM, but 
can be important if the corresponding Higgs couplings 
are enhanced. For Higgs production in association with 
top quarks, we find 
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Again, since the 2HDM Type MFV model shares the 
same modified top quark coupling as the 2HDM Types I- 
IV, they lead to a modified cross section that is invariant 
with respect to the Type of model. In our numerical 
results, we take from [51, 52]. 

Finally, we comment on the production cross section 
of the heavy CP-even Higgs H and the CP-odd Higgs A. 



For large couplings to bottom quarks, the main produc- 
tion mode is in association with bottom quarks. Here the 
cross section is in general completely different compared 
to 2HDM Types LIV, 



^Mi^ V ^BM „ i(tl,/\\ 

"bbH,A — "bbH,A X 1,?^ ) 



(21) 



In our numerical results, we use bbhOnnlo [53] to com- 
pute the SM cross section af^. The theoretical uncer- 
tainties of all SM production cross sections are also taken 
from [51, 52]. 

In summary, much of the light Higgs boson production 
cross sections in the 2HDM Type MFV are unchanged 
with respect to a 2HDM Type I-IV, with the notable 
exception of the production of the heavy Higgs bosons in 
association with bottom quarks. 



B. Decay Rates and Branching Ratios 

The partial widths of the Higgs bosons into fermions 
and weak gauge bosons can be written as 



YnffC^T^^fXitf? , (22) 

Tuvv ^ y^^v X {S.'irf , (23) 

where are the corresponding decay width of the SM 
Higgs boson. In our numerical analysis we compute these 
SM decay widths using HDECAY [54]. For the decays of 
the Higgs into gluons and photons we define the effective 
couplings 
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and compute at leading order 
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(24) 
(25) 



(26) 



As we only compute the ratio of partial widths, higher 
order corrections are expected to be small. To obtain ab- 
solute values of the decay widths into gluons and photons 
we use rf'^ from HDECAY [54] in Eqs. (24) and (25). 

Some of the lightest Higgs partial widths are Type- 
invariant, while others are Type-dependent. The Type- 
invariant widths include^ 



pMFV 

_ 1 1 

T^2HDM ^ ' r2HDM — ' 
^VV ^ gg 
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(27) 



At leading order, the width into diphotons depends on 
the Higgs coupling to the W boson as well as the top 
quark, which are both Type-independent couplings. In 



^ As already explained for the gluon fusion production cross section 
in Eq. (17), contributions from bottom loops can be neglected to 
a first approximation. 
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addition, all Types of 2HDMs can have a charged Higgs 
contribution to the diphoton and decay amplitudes. 
The size of the charged Higgs contribution depends on 
the the scalar trilinear coupling XhH±H± the charged 
Higgs mass. If the coupling and charged Higgs mass are 
the same between two different Types of 2HDMs, then 
the contribution is the same, giving 



- 1 , 



(28) 
(29) 



We must emphasize that the ratios, Eq. (29), are equal 
only when the charged Higgs coupling and mass are taken 
to be the same. As we will see, the range of the parame- 
ter space {^hH±H±7 Mff±) that is allowed by experimen- 
tal constraints may be considerably wider in the 2HDM 
Type MFV model versus the 2HDM Types I-IV, thus 
permitting larger effects on these rates. 

Finally, Higgs partial widths into fermions are Type- 
dependent, given by 
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C. Differences between the 2HDM Type MFV and 
Types I-IV 

In summary, there are two central differences between 
2HDM Type MFV versus Types I-IV: 

• The width r?|^^ can be completely different from 



bb 

(and r^M). 



p2HDM Ur^A rSM 



bb " bb 

• The width can be completely different from 

(and _). 

Since the total width of the light Higgs boson, Ftot, is 
dominated by F^j for = 125 GeV, this leads to the 
other important difference: 

• The total width F^f ^ can be significantly smaller 
or larger than F^hdm ^^^^ j^SMy 

The total width affects all of the branching fractions 
in a correlated way. To a very good approximation, 
Tr+T- ^ ^bb remains true for any 2HDM modification 
that can fit the ATLAS and CMS data. This is simply 
because the convolution of the production and decay of 
the Higgs boson to r+r" must be less than roughly the 
standard model rate. Hence, the total width ratio can be 
estimated as 



pMFV 
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The 2HDM ratios C^™'^ can be obtained by applying 
Eq. (10) to Eq. (9) for the couplings of the light Higgs h: 



IV. ONE LIGHT HIGGS BOSON 



d™ = ef™ = cosa/sin/3 (Type I) 
ef°^ = e.™ = -sina/cos/? (Type II) 



ef™ = cosa/sin/3, C. 



2HDM 



= — sin a/ cos f3 (Type HI) 



d = - sin a/ cos /3, = cos a/ sin /? (Type IV) . 

Here we see that the Type MFV interpolates among all 
other Types of 2HDMs following the limits in Eq. (10). 
There are, however, several fascinating parameter ranges 
that are not reached in any of these models. These in- 
clude: 



We now confront the 2HDM Type MFV with the avail- 
able data on standard model Higgs searches. We first 
consider a scenario with one light scalar Higgs boson at 
Mfi ~ 125 GeV. The second scalar boson H is assumed 
to be heavier than Mh > 150 GeV, such that it does not 
directly influence the interpretation of the data. A sce- 
nario with two light scalar bosons is discussed in Sec. V 
below. 



Fit to the Data 
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(32) 



Note that close to the limit ed,i — > — 1/ tan/3, the Yukawa 
couplings can become non-perturbatively large. 

Given that F(,{, is the dominant part of the total width 
of the lightest Higgs boson, these effects can have dra- 
matic consequences on all of the resulting Higgs branch- 
ing ratios. 



To determine the preferred values of the couplings of 
the 2HDM Type MFV given the current data, we perform 
a simple fit of Cvj Ct^i i and taking into account 
SM Higgs searches at LHC [1,2, 55] and at Tevatron [37]. 
Similar fits have also been performed in [28, 56-67]. We 
consider: 

(a) searches for h ^ bb produced in association with 
a gauge boson from ATLAS [68], CMS [69], and 
Tevatron [37], 

(b) searches for h ^ bb produced in association with 
top quarks by CMS [70] and ATLAS [71], 
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FIG. 1. Best fit regions in tiie - i,t (left), - (Jl (middle), and - (right) planes in a fit of the data to one Higgs 
boson at 125 GeV. The light green, green and dark green regions correspond to the A^^ =1, 4, and 9 regions. The red labeled 
contours in the left plot show constant values of tan /3. In the middle plot the region shaded in orange shows the parameter 
space that is accessible in a 2IIDM Type II by varying (^y within its Icr range. The blue solid curves in the middle and right 
plot exemplarily show regions of parameter space that can be reached in the MFV 2IIDM by varying within the la range, 
while keeping the other parameters fixed to the indicated values (dotted lines correspond to ^y outside the la range). For the 
blue curves in the right plot we fix tan/? = 0.78 and = —1.52 to the best fit values. The gray region in the left plot with 
<Q corresponds to a second minimum in the x^i that is however excluded by searches for the heavy scalar H (see text). 



(c) the inclusive h — > r+r searches at ATLAS [72] 
and CMS [73], 

(d) the inclusive h WW searches from ATLAS [74], 
CMS [75], and Tevatron [37], 

(e) the results on /i — )■ WW produced in VBF by 
CMS [76] and ATLAS [55], 

(f) the inclusive h ZZ — > 4£ searches from AT- 
LAS [77] and CMS [78], 

(g) the inclusive — > 77 searches from ATLAS [79], 
CMS [80], and Tevatron [37], and finally 

(h) the results on /i 77 produced in VBF by 
CMS [80] and ATLAS [55]. 

In a 2HDM the couplings of the Higgs boson to gauge 
bosons is constrained to be < 1. Furthermore certain 
regions of the parameter space are only accessible for very 
small values of tan /? that lead to non-perturbative and 
^~ couplings (see Eq. (9)). We therefore perform the fit 
imposing the constraints £,y < ^ and tan/3 > 0.5.'^ The 
resulting best fit regions of parameter space are shown in 
Fig. 1 in the - plane (left plot), the Sj^ - plane 
(center plot), and the ~ plane (right plot). The 
dark green, green and light green regions correspond to 
Ax^ =1,4, and 9, respectively, and we will refer to them 
as 1, 2 and 3(j regions. Throughout each plot the other 
couplings are chosen to minimize the total x^- 



We checked that none of our conclusion changes by allowing for 
even lower values of tan /3 > 0.3. 



In the best fit region, the coupling is to a good ap- 
proximation SM-like while the coupling is reduced but 
still positive. The reduced gives a slight enhancement 
of the partial width of — 77. Simultaneously, it also 
reduces the gluon fusion production cross section. There- 
fore, in order to obtain an enhanced inclusive h — >■ 77 rate 
as hinted by the data, a reduction of the total width is 
also required. This is achieved by reducing the cou- 
pling as shown in the center plot of Fig. 1. Both signs of 
the coupling are allowed and give essentially equiva- 
lent results for the light Higgs boson. Finally as shown in 
the right plot of Fig. 1, a reduced coupling is preferred 
because it leads to a strongly reduced h — > r+r^ signal, 
as hinted by CMS data [73]. 

The best fit values that we find read 

e^=0.99, Ct'* = 0.79, 
4^ = ±0.73, Cr=0. (34) 

We remark that there exists another minimum in the 
in a region where the coupling has the opposite sign 
compared with the SM. This region is shown in gray 
in the left plot of Fig. 1. In this region one could ex- 
pect a considerable enhancement of the /i — >■ 77 partial 
width as W and top loops now interfere constructively. 
However, since the Higgs coupling with W^s is generically 
suppressed, the enhancement of the di-photon width is 
small or even absent. Furthermore, as is evident from 
the left plot of Fig. 1, the opposite sign solution for 
requires tan /3 as small as possible. The requirement 
tan P > 0.5 then necessarily implies that the £,y coupling 
is considerably reduced in magnitude compared to the 
SM values. For tan/3 < 1 the sum rule, Eq. (11), im- 
plies that the heavy scalar bosons H has an enhanced 
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FIG. 2. Results for various Higgs rates normalized to the SM 
rates from a fit of the data to the 2HDM Type MFV with one 
light scalar boson at 125 GeV. For comparison results from 
an analogous fit in the 2HDM Type II and the experimental 
la ranges are also shown. The black stars correspond to the 
best fit values. The circles indicate an example scenario with 
strongly enhanced VBF /i — >■ 77 signal. 



gluon fusion production cross section. Furthermore, due 
to the sum rule, Eq. (14), both scalar bosons couple in a 
non-negligible way to vector bosons. Correspondingly, H 
has a large rate in H ^ WW/ZZ. Therefore SM Higgs 
searches rule out the region with < for masses up 
to Mh < 600 GeV, the present limit of the experimen- 
tal search sensitivities. For Mh > 600 GeV, the model 
is already in the decoupling limit, where we expects all 
light Higgs boson couplings to be SM-like. We will not 
consider this region further in this paper. 

Fig. 2 shows the rates of the light Higgs boson in the 
best fit region with > in comparison to the exper- 
imental data. The dark green, green, and light green 
bands correspond to the 1, 2, and 3cr regions of the fit. 
The black stars mark the best fit values. The rates ob- 
tained for the best fit values follow the data closely: the 
h — >■ T+T~ rate is reduced to near zero; the largest en- 
hancement is realized for /i — >■ 77 in VBF, followed by 
the inclusive h — )■ 77; the h — >■ WW and h — ?► ZZ rates 
are only slightly enhanced; the h ^ hb rate is reduced 



compared to the SM, which is in slight tension with the 
Tevatron data. Ignoring the Tevatron h ^ bh data in the 
fit would allow even larger /i -> 77 rates, by reducing 
h ^ bb further. 



B. Comparison with the 2HDM Type II 

In Fig. 2 we also compare the best fit values for the 
rates of the light Higgs in the 2HDM Type MFV with 
the corresponding rates in a 2HDM Type II. The most 
important differences are in the inclusive h — >■ t+t~ rate 
and the VBF ft — 77 rate. The r+r~ rate cannot be 
reduced to zero in the 2HDM of Type II in contrast to 
the 2HDM Type MFV. In a Type II model = Sj^ and 
therefore a strongly suppressed h — >■ r+r" rate implies a 
strongly reduced h ^ bb width. This in turn would lead 
to a drastic enhancement of all other branching ratios 
and correspondingly to /i — ZZ and h — )■ WW rates far 
above what is allowed by current data. In the 2HDM 
Type MFV, instead, and are independent from 
each other and is unconstrained by rates other than 

h — > T+T^. 

Concerning the VBF ft, — >■ 77 rate, we observe that 
in the MFV model huge enhancements are possible. In 
Fig. 2 we show an example scenario (open circles) lead- 
ing to an enhancement by a factor of ~ 3.5 that can be 
realized by the couplings 



0.97 



0.49 



^ = 0.33 



(35) 



By contrast, in the Type II model, we find the VBF 
ft — > 77 rate is bounded by approximately 1.7 at the 3a 
level. This difference can be traced back to the strong 
correlation between and in the Type II model. An 
enhancement of the VBF h —i' j-f rate by factors of a few 
is only possible if the h ^ bb width is reduced consider- 
ably, by reducing f^'. This not only enhances the ft — ?> 77 
branching ratio but also all other branching ratios. To 
keep the inclusive ft — > WW, ft — > ZZ, and ft — > 77 rates 
at a level compatible with experimental data, the gluon 
fusion production cross section has to be reduced by re- 
ducing the ^t'^ coupling. In contrast to the MFV model 
where and are independent, this is not possible in 
a Type II model, where a strongly modified coupling 
implies ^'^ ~ 1 and vice versa. 

We remark that an enhancement of the VBF ft — > 77 
rate by a factor of few in the MFV model also im- 
plies a similar enhancement of the VBF ft — )• WW rate. 
This is in tension with the CMS and ATLAS analyses 
of ft — WW —T' ivlv that do not see any excess above 
background in the ft + 2 jets sample that is dominated 
by VBF production. Note that these results are included 
in our fit, but given the considerable uncertainties, they 
do not influence the fit by much. Updated results for 
the VBF ft ^ WW rate from ATLAS and CMS can ei- 
ther rule out or support the possibility of an strongly 
enhanced VBF ft WW rate in the 2HDM Type MFV. 
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We also note that a strongly enhanced Wh — WWW 
rate is predicted in the regions of parameter space with a 
strongly enhanced VBF ft, — > 77 rate. However, the cor- 
responding CMS and ATLAS searches [81, 82] currently 
give only very mild constraints on this channel. 



C. Generic Predictions for the Light Higgs 







assoc. 
rh^bb 


incl. 
H-h^VV 




=> 


< 1 


> 0.9 


^^^h^JJ > 1.5 


=> 


< 0.8 


> 1.2 








assoc. 
rh^bb 




VBF . 




< 0.9 


> 0.3 


VBF ^ 




< 0.4 


> 0.7 



TABLE L Correlations between several Higgs rates and cross 
sections in the 2HDM Type MFV with one light Higgs boson 
at 125 GeV. 

In the 2HDM Type MFV, even though the couplings of 
the light Higgs to the different Types of fermions can be 
modified independently, there exists correlations among 
the rates of some of the currently investigated Higgs 
search channels. In Table I we summarize generic pre- 
dictions for the light Higgs that we found. 

In particular, we find that an enhancement of the in- 
clusive h 77 rate implies 

(i) an upper bound on the h ^ bb rate where the Higgs 
is produced in association with a vector boson, and 

(ii) a lower bound on the inclusive h -> WW and h — 
ZZ rates. 

Similarly, also strong enhancements of the VBF /i — 77 
rate imply stringent upper bounds on the h ^ bb rate. 
This is clearly shown in Fig. 2 where the circles corre- 
spond to an example scenario with the VBF h — > 77 
rate enhanced by a factor of ~ 3.5. As a very strong 
enhancement of the VBF /i -> 77 rate is only viable 
for a reduced coupling, the sum rule, Eq. (11), im- 
plies a lower bound on the production cross section of 
the heavy scalar H, and correspondingly good prospects 
for H searches at the LHC. Note that the bounds pre- 
sented in Table I do not change appreciably if we restrict 
to scenarios fitting the Higgs data at the 1, 2 or 3cr level. 
Correspondingly, these bounds are robust. 



D. The Quasi-Decoupling Limit 

We analyze the extent to which the identified values 
of the Higgs couplings can be realized concretely in the 
2HDM Type MFV. In the best fit region with positive, 
the coupling of the Higgs with gauge bosons is approx- 
imately SM-like. We can express this "quasi-decoupling 
limit" with the relation a = /3 — Tr/2 -\- x where a; is a 
small expansion parameter. The couplings of the Higgs 
bosons with gauge bosons and fermions are then given 

by 

CO=^l-y, (36) 
^v^x, (37) 

e::^ (^l-y)-fxC;f , (38) 
(l-y)+^, (39) 
e,^^(^l-y)-xe,^, , (40) 

C,^.^ (l-y)^.^. + ^, (41) 

with the couplings given in Eq. (9) . We see that even if 
the coupling of h to weak gauge bosons is SM-like within 
a few percent, there can be substantial modifications to 
the remaining couplings (this has been also pointed out 
recently in [29] in the context of the 2HDM Type I.). In 
particular we observe that couplings of the light Higgs 
with quarks and leptons can be strongly modified if the 
couplings of the CP-odd Higgs boson with quarks and 
leptons are sizable. In this situation the couplings of the 
heavy CP-even Higgs boson H are to a good approxima- 
tion the same as the CP-odd Higgs boson A.^ We again 
stress that in the MFV framework, the couplings of h to 
up-type quarks, down-type quarks, and leptons can be 
modified independently, whereas in 2HDM Types I-IV 
they are strongly correlated. For example, in the 2HDM 
Type II, an enhancement of = S,^ — tan (3 is always 
accompanied by a reduction of — l/tan/3. 

It is impressive that, in contrast to 2HDM Types I-IV, 
the 2HDM Type MFV can exactly reproduce the best fit 
values for the Higgs fermion couplings while retaining es- 
sentially SM-like couplings of the Higgs to gauge bosons. 
This is illustrated in the middle and right plot of Fig. 1. 
The additional parameter eb in the MFV framework al- 
lows us to modify and completely independently. 
Example choices for tan/3 and ef, that cover the whole 
parameter space of top and bottom couplings, and in 
particular to reach the best fit values, are indicated with 



^ Conversely, if the couplings to quarks or leptons of A are small, 
then the corresponding couplings of H can be modified, with the 
couplings of h remaining SM-like. 
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the blue solid lines in the middle plot. Similarly, the right 
plot of Fig. 1 shows how the independent parameter 
can be used to obtain a highly suppressed coupling 
while fixing tan /3 and e?, such that i^^'' and correspond 
to their best fit values. 

For completeness we also report our best fit values for 
the 2HDM parameters tan/3, a, ct, and e^.. We find 

tan/3 0.78 , a ^ -1.05 , 
eb = -1.52 (-8.3) , er = -1.74 , (42) 

where the two numbers for e^, correspond to the negative 
(positive) solution for In the best fit point, tan/3 
has a small value. From the left plot of Fig. 1 we note 
however, that also regions of parameter space with large 
tan /3 can result in a very good fit of the data. 

E. Predictions for the Heavy Scalar H 

Once the couplings of the light Higgs boson are fixed, 
then the couplings of the heavy Higgs boson H (as well as 
the CP-odd Higgs boson A and the charged Higgs bosons 
H^) are determined. At the best fit values, we find 

= 0.14 , = 1-36 , 
= -12.4 (-1.78) , = -7.1 . (43) 

In Fig. 3 we show the resulting predictions of various 
cross sections of the heavier scalar H as a function of 
its mass for the two solutions for ^^.^ The top (cen- 
ter) plot corresponds to the best fit point with > 
(^^ < 0). In the bottom plot we allow for a larger 
coupling such that the inclusive h — > r+r^ rate is 
50% of the corresponding SM rate^. Plotted are the cur- 
rent exclusion bounds normalized to the predicted signal 
cross sections. We take into account not only searches 
for SM-like Higgs in the H ^ ZZ ^ M [77, 78] (light 
blue), H ^ ZZ ~> [83, 84] (dark blue), and H 
WW -> dviv [75, 85] (green) channels, but also searches 
for MSSM Higgs bosons in the H t+t- [86, 87]^ (dark 
red) and H ^ bb [88, 89] (orange) channels. For every 
value of Mji we consider the strongest of the individual 
bounds from ATLAS and CMS. 

We find that for the best fit point with > (top 
panel), the heavy scalar is excluded up to AIh ^ 350 GeV 
by the current searches. Once Mh > 350 GeV, the heavy 
scalar has a large branching fraction into tt and is cor- 
respondingly only weakly constrained by present data. 



tan/? = 0.78 , = -1.05 , e,, = -8.3 , e, = -1.74 
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FIG. 3. Current 95% C.L. exclusion bounds normalized to 
the predicted cross sections for the heavy Higgs as function of 
the heavy Higgs mass Mh- Shown are H — ^ ZZ — >■ 4^ (light 
blue), H ^ ZZ ^ llvu (dark blue), H WW ivtv 
(green), H — >■ t^t" (dark red), and H ^ bb (orange). The 
top (center) plot corresponds to the best fit values of the 
light Higgs coupHngs with > (^^ < 0). In the bottom 
plot we allow for a larger coupling such that the inclusive 
h — >■ T^T^ rate is 50% of the corresponding SM rate. 



^ Here we do not consider the possibility of a sizable H hh decay 
rate, as the corresponding coupling is another free parameter of 
the model. 

® As shown in Fig. 2, — > r+r" rates around 50% are still result- 
ing in a reasonable good fit of the present Higgs data. 
Note that the CMS analysis [86] only provides bounds in the 
M^-tan /3 plane of a specific MSSM scenario. We translate these 
bounds into bounds on the signal cross sections and reinterpret 
them in our scenario assuming constant efficiencies. 



In the < case, the coupling of the heavy Higgs to 
bottom quarks is considerably larger, c.f. Eq. (43). As a 
result, the branching ratios of H into t^t^ and vector 
bosons are significantly smaller and current searches only 
start to be sensitive to the region below ^ 350 GeV. 
Searches in the di-tau final state are able to just barely 



9 



exclude the heavy scalar up to Mh — 200 GeV. Searches 
for h — ZZ — )■ M start to become sensitive to H in the 
mass range up to 350 GeV. Above 350 GeV, the heavy 
scalar again decays dominantly into tt. We remark that 
the reason for the strong constraints from searches in the 
T+T^ final state is the best fit preference to suppress 
the h T^T^ rate as much as possible. Indeed a 
coupling close to zero is only possible with an enhanced 
coupling of the heavy scalar H to r's [see Eq. (13)]. If 
we deviate slightly from the best fit point and allow a 
somewhat larger h — )• r+r^ rate, the parameter space 
for the heavy scalar opens up further. This is shown in 
the bottom panel of Fig. 3 where we fix ei such that the 
inclusive h — > r+r^ rate is 50% of the corresponding 
SM rate. Now the coupling is much smaller and the 
heavy Higgs is only excluded in a small region around 
300 GeV by searches in the Ai final state. Interestingly, 
apart from the h — > t~^t~ rate, all other channels are 
to an excellent approximation unaffected by changing ei. 
Note also that in this scenario the current H ^ bb con- 
straints are stronger than the H — >■ t+t^ ones, show- 
ing that H ^ bb MSSM searches at the LHC [90] give 
valuable complementary information in the 2IIDM Type 
MFV. 

If the other couplings of the light Higgs are modified 
from their best fit values, prospects to probe the heavy 
Higgs typically remain excellent throughout large parts 
of the parameter space: 

• For a closer to 1 and correspondingly for larger 
values of tan /3, the gluon fusion production cross 
section for H can be reduced. Still even for large 
tan /3 we find that the heavy Higgs signals are typ- 
ically only a factor of a few above the current ex- 
clusion limits. 

• For larger tan f3 the coupling of H to the top quark 
is reduced, implying reduced branching ratios of 
H ^ tt for Mh > 350 GeV and therefore in- 
creased sensitivity of iJ — >■ WW/ZZ in the large 
Mh regime. 

• A reduced coupling leads to a larger coupling 
and therefore the H — WW and H — )• ZZ signals 
get enhanced. The current Higgs searches already 
probe the corresponding parts of parameter space. 

• As already discussed, the value of the coupling 
strongly influences the branching ratios of H. The 
larger the deviation of from 1, the larger is the 

coupling [see Eqs. (40), (41) and the discussion 
below]. Nonetheless, we find that even the largest 
deviations of from 1 generically lead to signals 
in the H WW/ZZ channels that can observed 
in the near future, as long as Mh ^ 350 GeV. 

There also exist also corners of parameter space where 
the heavy scalar cannot be detected. If the couplings of 
the light Higgs to the top quark and gauge bosons are 
to a very high precision SM-like and tan /3 is large, the 



couplings of H to the top quark and gauge bosons can be 
made arbitrarily small. This results in very small produc- 
tion cross sections for H that are not easily detectable. 

The phenomenology of the heavy CP-odd Higgs A dif- 
fers from that of the CP-even H, as A does not couple 
to weak gauge bosons. Thus the only search channels 
that are currently able to probe A are the A — >■ r+r" 
and A —i' bb searches. Analogous to the H boson, for 
Ma 350 GeV, the A —i' tt decay can open up. 

V. TWO LIGHT HIGGS BOSONS 

We now investigate a scenario where, in addition to 
the lightest Higgs boson at 125 GeV, the second CP even 
Higgs H boson is also light. For definiteness we fix its 
mass to Mh = 135 GeV. We checked that varying the 
mass of the second Higgs between 130 GeV and 140 GeV 
does not change the results qualitatively. For a recent 
discussion of a similar scenario in the context of the 
NMSSM see [36]. Scenarios with two quasi-degenerate 
Higgs bosons at 125 GeV are discussed in [33-35]. 

A. Fit to the Data 

In the Higgs search channels with mass resolution 
smaller than Mh — Mh — 10 GeV, specifically h — > bb, 
h ^ T^T~ and h — WW, the LHC collaborations would 
effectively be observing the sum over the signals coming 
from both Higgs bosons h and H. For the high resolution 
channels /i — 77 and h — > ZZ, we instead consider the 
data on the signal strength separately for h and H. 

The resulting best fit regions in the - - 
and - planes are shown in Fig. 4. Interestingly 
enough, now that the effects of the second Higgs boson 
are directly included in the fit, in the £,y - plane only 
one solution with ^j/' > is present. Both couplings are 
slightly reduced compared to the SM values. The two 
signs of the coupling are not equivalent anymore. We 
find a slight preference for a negative that is below 1 
in magnitude. The preferred coupling remaining close 
to zero. The best fit point is given by 

eO=0.85, Ct'*-0.77, 
^fc = -0.52 , = 0.16 . (44) 
The resulting couplings of the heavier Higgs are 
e^ = 0.53, ef = 0.66, 
ef--2.7, ef--1.6. (45) 

The reduced couplings of the heavier Higgs to the 
top quark and vector bosons lead to suppressed pro- 
duction cross sections of H. The enhanced coupling to 
bottom quarks, however, results in a branching ratio of 
BR{H — >• bb) ~ 95%. Correspondingly, the second Higgs 
primarily just adds to the h/H — ?> bb signal while its ef- 
fect in the other search channels is negligible to a first 
approximation. 
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FIG. 4. Best fit regions in the (,v - Ct" (left), Cb " it (middle), and - (right) planes in a fit of the data to two Higgs 
bosons at 125 GeV and 135 GeV, respectively. The dark green, green and light green regions correspond to the Ax^ =1, 4, 
and 9 regions. The red labeled contours in the left plot show constant values of tan /?. In the middle plot the region shaded 
in orange shows the parameter space that is accessible in a 2HDM Type II by varying (^y within its la range. The blue solid 
curves in the middle and right plot exemplarily show regions of parameter space that can be reached in the MFV 2IIDM by 
varying within the la range, while keeping the other parameters fixed to the indicated values. For the blue curves in the 
right plot we fix tan/3 = 6.7 and eb = —0.57 to the best fit values. 



The best fit values can be accommodated in the 2HDM 
Type MFV for appropriate choices of Cf, and eg. Example 
choices for tan e;,, and eg that allow to cover the whole 
parameter space of top, bottom and tau couplings, and 
in particular that reach the best fit values, are indicated 
with the blue solid lines in the middle and right plot of 
Fig. 4. In the right plot, tan /3 and are fixed such that 
and correspond to their best fit values. We obtain 
the best fit values 

tan/3 = 6.7 , a = -0.71 , 

ef, = -0.57 , Er = -1.0 . (46) 

By contrast, in the 2HDM Types I-IV, the best fit val- 
ues cannot be accommodated, even if these models can 
still produce a reasonable good fit of the present LHC 
Higgs data. The two bands shaded in orange in the mid- 
dle plot of Fig. 4 show the region that can be reached 
in a 2HDM Type II, if the coupling is varied in the 
1(7 range around the best fit value. This region does not 
cover the best fit values for and 

The Higgs rates resulting from the fit of the two light 
Higgs boson scenario are compared to the experimental 
data and to the corresponding fit in the 2HDM Type II 
in Fig. 5. The best fit points are shown as black stars, 
while the 1, 2, and 3cr ranges are indicated by the green 
shaded bands. In the best fit point the /i — )■ 77 and 
h — VV rates are slightly enhanced due to the suppres- 
sion of the coupling that controls the total width of 
h. The dramatic reduction of the coupling leads to 
a /i — > T+T^ rate close to zero. Note however that in 
contrast to the one Higgs case, the r+T~ rate cannot 
be eliminated completely, because the second scalar H 
always contributes at some level. The main difference 
from the single light Higgs boson case (Sec. IV) is that 



the /i — > 66 rate is slightly enhanced compared to the SM 
prediction, with h and H decays contributing approxi- 
mately 65% and 45% of the SM rate. At the ia level, 
the bb rate can even be enhanced by up to a factor of 
1.6. Overall, we find that the fit with the two light Higgs 
bosons is just as good as the scenario with only one light 
Higgs boson. 



B. Generic Predictions for Higgs Signals 
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TABLE II. Correlations between several Higgs rates in the 
2IIDM Type MFV with two light Higgs bosons at Mh = 
125 GeV and Mh = 135 GeV. 

We now discuss generic correlations for the Higgs sig- 
nals in the 2HDM Type MFV with two light scalar 
bosons. As shown in Table II, we find that just like 
the one light Higgs boson scenario, an enhanced inclu- 
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FIG. 6. Values of the /i — >■ 77 rate, normalized to the SM 
rate in the Mfj± - \hH±H± plane. Tree level couplings of the 
light Higgs to gauge bosons and fermions are assumed to be 
SM-hke (C^ = = !)• 



VI. IMPACT OF THE CHARGED HIGGS 
BOSON 



FIG. 5. Results for various Higgs rates normalized to the 
SM rates from a fit of the data to the 2HDM Type MFV 
with two light scalar bosons at 125 GeV and 135 GeV. For 
comparison, results from an analogous fit in the 2HDM Type 
II and the experimental la ranges are also shown. The black 
stars correspond to the best fit values. 



sive /i — ?> 77 rate implies: 

(i) an upper bound on the h + H ^ bb rate where the 
Higgses are produced in association with a vector 
boson, and 

(ii) a lower bound on the inclusive h + H ^ WW rate. 

Quantitatively however, because of the additional con- 
tribution of the second scalar H, both the upper bound 
on the bb rate and the lower bound on the WW rate are 
slightly larger than in the one Higgs case. Analogously, 
the upper bounds on the h + H ^ bb rate that are im- 
plied by strong enhancements of the VBF ft, — > 77 rate 
are weaker compared to the one Higgs case. 

Note that in this setup a second peak at the mass of 
the heavier Higgs is expected in the — > 77 and H — >• 
ZZ channels. However, given the suppressed and 
couplings, the corresponding signal strengths are often 
only few percent of the SM signals. Therefore, finding 
evidence for a second light Higgs in these channels is 
challenging. 



A well-known additional non-standard effect in the 
phenomenology of the light Higgs boson can arise from 
its interaction with a light charged Higgs that can con- 
tribute at I-loop to the decay /i — ?> 77 (for studies of loop 
induced corrections to the /i — > 77 rate in new physics 
models see e.g. [27, 56, 91-103]. The main contributions 
to the partial width of the Higgs into two photons are 
then coming from W loops, top loops, and charged Higgs 
loops 



T{h -> 77) 



2567r3 i>2 



(47) 



where x, = AMf/Ml, ^ 3, Qt ^ 2/3, and \hH±H± is 
the coupling of the light Higgs boson with two charged 
Higgs bosons divided by v. Finally Ai, A1/2 and are 
loop functions of a gauge boson, a fermion, and a scalar, 
respectively, and are collected in Appendix B. 

As shown in Fig. 6, sizable corrections to the h — > 
77 partial width are only possible for very light charged 
Higgs bosons and for large negative coupling with the 
lightest Higgs boson. 

First consider the familiar case of the MSSM, where 
the coupling of the lightest Higgs boson with the charged 
Higgs is determined by electroweak gauge couplings and 
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erf=0,e, = Ed = -1-52, = -1.74 = -0.57 , e, = -1. 




(GeV) Mfi' (GeV) M„. (GeV) 



FIG. 7. Constraints in the M^fi - tan/3 plane from direct searches for charged Higgs bosons in top decays (the orange 
region/solid contour corresponds to the H ^ ti/ final state, the yellow region/dotted contour corresponds to the H ^ ca final 
state) and from measurements of the B tv branching ratio (red region/dashed contour). All constraints are 95% C.L. 



is given by 



Ai -|- A2 



UMSSM I 
\'^hH±H± 



9I +9I .9I 



< 52 

^ 2 



0.21. 



(48) 



This coupling is too small to give a visible effect on the 
h — > 77 partial width. 

On the contrary, in a generic 2HDM there exist regions 
of parameter space producing a large negative coupling 
^hH±H± ~ —3 and light charged Higgs bosons, while 
still remaining compatible with constraints from vacuum 
stability and electroweak precision observables. These 
regions can produce an enhancement of the /i — > 77 rate 
by a factor as large as 1.5-2. 

The full expression for Xhh^h^ ii^ terms of a gen- 
eral 2HDM Higgs potential parameters is given in Ap- 
pendix A. However, it is easier to understand the param- 
eter dependence of XhH±H± and of the charged Higgs 
mass in the "almost" decoupling limit and at large tan (3. 
In this regime, using the definition of the scalar potential, 
Eq. (Al), we have 



XhH±H± ^ (A3 + A4) . 



(49) 



where it follows that A3, A4 must be negative and sizable 
to have a sizable enhancement of the Higgs to di-photon 
rate. Additionally, A4 < allows a smaller mass for the 
charged Higgs boson given that 



Mi + (A4 + A5)- 



(50) 



and thus further increases the /i — > 77 rate. 

Large negative couplings in the general 2HDM poten- 
tial could lead to vacuum stabilities issues. Following 
[104], vacuum stability at tree-level is satisfied so long as 

Al, A2 > 0, 



+ A3 + A5-2IA6 + A7I >0, 



A3 > -V A1A2 , A3 — IA5I > — ■\/AiA2 



(51) 



in order to have a potential bounded from below. A full 
treatment of vacuum stability requires minimization of 
the renormalization-group improved potential, but this 
is beyond the scope of this paper. Large negative A.^ cou- 
plings could lead to additional minima deeper than the 
electroweak breaking minimum, but we checked at tree- 
level that this does not occur for the scenarios considered 
here. Finally, we have also checked that our Higgs po- 
tential parameters satisfy electroweak precision observ- 
ables. We find that Ai* < 0.1 throughout the parameter 
space, while AT < 0.1 so long as either Mh± ~ Ma or 
Mh± ^ Mh (or both). We find that a large coupling 
XhH±H± ^ ~3 is compatible with all constraints, as long 
as A5 7^ 0. 

Very light charged Higgs bosons are subject to several 
constraints. Model independent bounds from LEP ex- 
clude charged Higgs masses below < 80 GeV [105]. At 
the Tevatron and the LHC, charged Higgs bosons are 
searched for in decays of top quarks, with the charged 
Higgs decaying either to tv or into two jets. Limits are 
obtained for the product of the branching ratios BR(t — ?> 
i?±fo)xBR(i/± ^ Tv) and BR{t i?±fe)xBR(i/± ^ 
jj). Additional constraints arise from flavor observ- 
ables, in particular from the branching ratio of the decay 
Bu TV that is sensitive to the tree-level exchange of a 
virtual charged Higgs. In a 2HDM Type II, all these lim- 
its are conveniently presented in the A/^± - tan/3 plane: 
the only two parameters that the branching ratios depend 
on. 

The left plot of Fig. 7 shows a summary of these con- 
straints in the Type II model. The orange solid contour 
corresponds to a combination of the t — >■ hH^ — ?> tv 
searches at ATLAS [106], CMS [107], and DO [108]. 
The yellow region inside the dotted contour is excluded 
by a combination of t — ^ hH^ — >■ cs searches at AT- 
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LAS [109], CDF [110], and DO [108]. The red region 
above the dashed contour is excluded by the latest com- 
bination of Bu — ?> Ti/ data from BaBar [111, 112] and 
Belle [113, 114], that is in reasonable agreement with SM 
expectations. The bounds from top decays only exist for 
AIh± 160 GeV, while obviously no such restriction ex- 
ist for the constraint from i?„ — tv. We observe that 
in the Type II model a window around tan /3 ~ 6 and 
Mfj± ~ 100 GeV cannot be excluded based on current 
available data. In this region of parameter space, charged 
Higgs loops can lead to large enhancements of the /i — 77 
rate. 

In the 2HDM Type MFV, the constraints in the Mh± ~ 
tan /3 plane depend strongly on the e parameters. Gener- 
ically we find that, large values of tan/3 > 5 are mostly 
unconstrained by current data for factors of 0(1). Cor- 
respondingly, in the 2IIDM Type MFV, large regions of 
parameter space are open where charged Higgs loops can 
enhance the /i — >■ 77 rate significantly. This is illustrated 
in the center and right plot of Fig. 7, that show again the 
constraints in the Mj:/±-tan/3 plane, fixing the values for 
£(, and Cr to the best fit values in the one light Higgs case 
(center) and the two light Higgs case (right) as indicated. 
Note that these values represent simply example scenar- 
ios, since, introducing the effects of the charged Higgs in 
the di-photon rate, the best fit values will change. Low 
values of tan /? generically remain constrained by charged 
Higgs searches in top decays. Indeed, low values of tan/3 
correspond to sizable tj^bLH^ couplings and therefore to 
large production from top decay. 

In principle also the loop induced — > 57 decay sets 
strong constraints on light charged Higgs bosons. For 
example, in the 2HDM Type II, the bound Mh± > 
380 GeV holds [115, 116]. This bound would rule out 
visible charged Higgs effects in ft — )• 77. However, go- 
ing beyond the Type II model, the 6 — > S7 bound de- 
pends not only on the charged Higgs mass, but also on 
the charged Higgs couplings and [117, 118]. More- 
over, being a FCNC process, 6 57 is also sensitive to 
higher order terms in the expansions of the Higgs cou- 
plings, Eqs. (3), (4), that are not relevant for Higgs col- 
lider phenomenology. Correspondingly, we do not con- 
sider constraints from 6 ^ 57 here. 



VII. CONCLUSIONS 

We have explored the detailed Higgs phenomenol- 
ogy of a 2HDM based on the MFV principle, in which 
both Higgs doublets couple to up-type and down-type 
fermions. The agility of the model permits several pos- 
sibilities to explain the current hints of ft. — >■ 77 excess 
in both the inclusive and exclusive VBF channels, with- 
out leading to large enhancements in the gauge boson or 
fermionic channels. Simultaneously, the model allows to 
accommodate a strongly reduced ft — > t+t^ rate, still 
having a SM-like h bb rate. 

Current Higgs data is well described in two distinct 



regimes of the light Higgs couplings: One regime is where 
the coupling of the Higgs to top quarks has opposite sign 
with respect to the Higgs - gauge boson couplings, and 
the second is a quasi-decoupling regime. We find that the 
first regime is ruled out by searches for the heavy Higgs 
boson H. The quasi-decoupling regime is the only viable 
region of parameter space in the 2HDM Type MFV. 

In this regime we find the VBF ft — > 77 can be en- 
hanced by up to a factor of 3 or more above the SM 
rate, still being consistent with the present LHC Higgs 
data. This occurs by simultaneously: reducing Ff,^ sub- 
stantially through the MFV parameter e^; virtually elim- 
inating Tt-+t.- ~ through e^; reducing the gg ^ h 
production through a reduction in the coupling of the 
Higgs with top quarks = cos a/ sin /3; while leaving 
the coupling to gauge bosons nearly identical to the SM, 
~ 1. As we showed in Table I, a large enhancement in 
the VBF di-photon channel has important consequences 
in the phenomenology of the light as well as of the heavy 
Higgs. An enhanced VBF ft ^ 77 rate implies both a 
suppressed ft — &6 rate and a lower bound on the pro- 
duction cross section of the heavy Higgs boson H. In 
particular an enhancement in the VBF di-photon chan- 
nel by a factor of 2-3 would automatically imply very 
good prospects for the detection of the heavy Higgs bo- 
son. 

The inclusive ft — ^ 77 rate can be enhanced by up to a 
factor of ~ 1.5, purely through a suppression of the total 
light Higgs boson width. Again there are correlations be- 
tween an enhanced inclusive ft — > 77 rate with a slightly 
suppressed ft — )■ 66 rate and a slightly enhanced ft — VV 
rate, as shown in Table I. Here it is important to point 
out that the Type MFV model provides a better fit to 
the existing (combined) data compared to a 2HDM Type 
II, mainly because the prediction for ft — ^ VV rates can 
be lowered slightly relative to a Type II 2HDM. 

In addition to the width effects, we also showed that 
a light charged Higgs boson can lead to significant loop- 
induced enhancements in the decay rate ft — >■ 77, up to a 
factor of 2 relative to the SM rate, as shown in Fig. 6. For 
a given charged Higgs mass and coupling \hH±H±i the 
contribution is otherwise the same between the 2HDM 
Type MFV versus 2HDM Types I-IV. However, there is 
a considerably wider range of (Af^±, XhH^H^) that is 
permitted in the 2HDM Type MFV. When the charged 
Higgs is light enough for t — >■ if +6 to be present, much of 
the parameter space of a 2HDM Type II model is ruled 
out by constraints on this rare top decay mode. In addi- 
tion, B — > also rules out a large swath of parameter 
space at larger tan (3 > 5-20 for charged Higgs masses be- 
tween 80-200 GeV. The 2HDM Type MFV is far less re- 
stricted by these constraints. Generally, for ed.e = 0(1), 
the region above tan/3 > 3 is fully allowed. Moreover, the 
usual difficulties of accommodating such a light charged 
Higgs boson from the constraints on 6 57 can be miti- 
gated given that this flavor-changing process is also sen- 
sitive to higher order terms in the MFV expansion for 
the Yukawa couplings shown in Eqs. (3), (4), that were 
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not relevant for Higgs collider phenomenology. 

A quasi-decoupling regime provides the best fit the 
present 125 GeV Higgs data. The deviations from the 
exact decoupling limit, that seem to be hinted by the 
present data, imply excellent prospects for observing the 
heavier Higgs H in channels typically used to search for 
a heavy SM-like Higgs boson. 

This we showed in Fig. 3, where we find that, corre- 
sponding to scenarios that are good fits to the present 125 
GeV Higgs data, the existing SM-like searches are typi- 
cally sensitive to the heavier Higgs up to about Mu < 
350 GeV. In addition, the search for H ^ bb can in some 
cases provide stronger constraints than H — >• r+r", due 
to the enhanced coupling of H to 6-quarks. Above this 
value for the heavier Higgs mass, the mode H tt opens 
up, and can have a large branching fraction that is only 
weakly constrained by present data. This also provides 
an exciting opportunity - the possibility of searching for 
the heavier Higgs of a 2HDM Type MFV through a t-t 
resonance. 

The properties of the heavy Higgs boson typically differ 
substantially from the heavy Higgs bosons of the MSSM. 
The overwhelmingly dominant decay modes throughout 
(essentially) the full mass range are H bb and H — 
T+r~, where dedicated searches exist. We emphasize 
that in order to probe the heavy scalar of the 2HDM Type 
MFV, and in general to find or rule out general 2HDMs, 
it is essential to continue the Higgs searches across all 
search channels in particular in SM-Higgs channels in 
mass regions where a SM-like Higgs is already ruled out. 

Finally, we showed that both h and H could be si- 
multaneously light, which can provide qualitative differ- 
ences in the SM-like Higgs rates that have comparatively 
low mass resolution. This means that, when h and H 
are within ~ 10 GeV of one another, some rates can 
be added, including h + H ^ WW, h + H ^ bb, and 
h + H — > r+T~. This "two light Higgs boson" case means 
the total rate into bb could exceed the SM value, while 
simultaneously having a sizable increase in the inclusive 
and exclusive VBF rates of /i — )■ 77. This is generally 
not possible in the case of just one light Higgs boson. 
In addition, in this scenario, there is a small rate of H 
into the high resolution channels H — ?> 77 and H — >■ ZZ. 
The continued exploration for resonances in these SM- 
like Higgs search channels, but at suppressed rates, would 
provide an outstanding opportunity to find or rule out 
this intriguing two light Higgs boson framework. 
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Appendix A: Two Higgs Doublet Model Potential 

The most general renormalizable Higgs potential of a 
2HDM can be written as 



(Al) 



Ai 



A5 



h.c. 



with {H2H1) = H+H^ - H^H°. The parameters Bfi, 
A5, Ag and Ay are in general complex. However in our 
study we consider all the parameters real. 

Assuming CP conservation, the Higgs fields entering 
Eq. (Al) can be parameterized in the following way 



Ho 



75(^^/5 



.V2 



ia2) 



(A2) 



and the real and imaginary parts of the Higgs fields do 
not mix. The physical spectrum consists of a charged 
Higgs iJ*, two scalars h and H and a pseudoscalar A. 
The Goldstone bosons and G provide the longitudinal 
components of the W and Z bosons, respectively 




(G\ Is0 -cp\ (a2\ 
\AJ sp ) \ai) 

(h^) - (:: 7;) (S ■ 

The coupling of the light scalar h with two charged 
Higgs bosons can be written in terms of the A.; couplings, 
tan /3, and a as 

Kh±h± = -AiSqS^C/3 + \2Cacls13 (A4) 

+ ^eica+psl + 2spSacf^) + XriCa+pCp + 2cpCasf)) . 

Appendix B: One Loop Functions 

The loop functions appearing in the expression for the 
/i 77 partial width, Eq. (47), are given by 



Ai(x) 



[2x-'^ + 3x~^ + 3(2x~i - l)f{x-^)] 
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Ai/2(.t) = 2x2 [x-' + {x'^ - l)f{x-')] 
Ao{x) = -x^ [x-' - fix-')] , 



with f{z) = arcsin^(\/^) for z < 1, which is the case that 
is relevant for us. 
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